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I.  MODELING  OF  IONIZED  GAS  SYSTEMS 
Program  Leader  -  Dr.  A.  V.  Phelps 

Our  studies  of  the  modeling  of  Ionized  gas  systems  include  a)  the 
theoretical  prediction  o‘f  the  conditions  for  the  growth  of  a  constriction 
In  a  weakly  ionized  gas,  such  as  found  in  a  molecular  gas  laser;  b)  the 
development  of  theories  of  radiation  transport  and  hydrodynamics  applic¬ 
able  to  a  wider  variety  of  real  situations,  such  as  expanding  and  turbu¬ 
lent  plasmas  in  geometries  other  than  the  usual  planar  geometry;  c)  the 
testing  of  radiative  transport  theories  under  controlled  laboratory  con¬ 
ditions;  and  d)  the  experimental  determination  of  the  laws  of  radiative 
scattering  and  absorption  appropriate  to  radiative  transport  problems. 

The  overall  objective  of  this  program  is  to  supply  experimentally  tested 
tinders tanding  and  computational  techniques  required  for  the  analysis  and 
-prediction  of  the  role  of  a)  the  generation  and  transport  of  radiation 
in  plasmas  and  b)  the  transport  and  reactions  of  charged  particles  in  weak 
to  fully  ionized  gases.  Because  of  its  concern  with  a  number  of  simul¬ 
taneous  atomic  processes  and  with  relatively  large  numbers  of  interacting 
atoms  or  molecules,  this  program  also  provides  an  important  link  between 
-the  experimental  and  theoretical  determinations  of  cross  sections  and  rate 
coefficients  which  make  up  much  of  the  remainder  of  the  program  and  the 
technology  of  ionized  and  excited  gases  which  this  work  is  designed  to 
support. 


Sterility  of  Discharges  In  Weakly  Ionized  Gases  (Drs.  E.  F.  Jaeger, 

L.  Os ter,  and  A.  V.  Phelps) 

In  this  theoretical  study,  numerical  calculations  have  been  made 
the  transient  behavior  of  an  infinitely  long,  cylindrical  discharge  in 
helium  at  pressures  from  100  to  1000  Torr.  Such  calculations  are  of 
interest  in  understanding  the  stability  of  weakly  ionized  gas  discharges 
such  as  occur  in  high  power  molecular  lasers. 

In  contrast  with  previous  calculations  ,  full  account  is  taken  of 
the  specific  heat  of  the  neutral  gas,  electronic  recombination,  and  the 
variation  of  gas  pressure  in  both  space  and  time.  Starting  with  the 
number,  momentum,  and  energy  balances  for  electrons,  ions,  and  neutrals, 
we  reduce  the  problem  to  a  set  of  five  time  dependent  partial  differential 
equations.  These  include  the  equation  of  thermal  energy,  the  ambipolar 
diffusion  equation,  the  continuity  equation  for  neutrals,  and  the  flux 
equations  for  charged  and  neutral  particles.  We  assume  that  the  discharge 
Is  operated  in  an  external  circuit  with  a  constant  voltage  source  and 
impedance.  In  the  numerical  solution,  the  axial  electric  field  in  the 
discharge  is  recomputed  after  each  time  step  to  account  for  the  decreasing 
voltage.  The  ionization  and  recombination  rater;  along  with  the  transport 
coefficients  for  helium  have  been  taken  as  functions  of  the  axial  electric 
field  divided  by  the  gas  density  from  the  best  experimental  measurements 
currently  available. 

*G.  Ecker,  W.  Kroll,  K.  H.  Spatschek,  and  0.  Zoller,  Physics  Fluids  10, 


1037  (1967). 


Figs.  1  and  2  illi-strate  the  numerical  solution  for  helium  initially 
at  100  Torr  pressure  and  300°K  with  an  initial  axial  electric  field  of 
400  V/cm.  The  external  resistance  is  100  ohms  per  cm  length  of  discharge. 
The  solution  is  started  by  turning  on  the  electric  field  at  time  zero 
in  the  presence  of  an  initial  charge  density  distribution  which  is  a 
Bess-1  function  falling  from  10 10  cm  3  on  axis  to  zero  at  the  wall.  The 
solid  lines  in  Fig.  1  show  the  transient  behavior  of  the  charge  density, 

“*  the  neutral  density,  nn>  and  the  neutral  temperature,  Tn,  computed 
at  the  center  of  the  discharge  for  both  10  and  50  grid  points.  We  find 
that  a  spatially  uniform  pressure  approximation  (dashed  lines  in  Fig.  1) 
gives  results  not  far  different  from  the  complete  solution.  In  Fig.  2, 
we  compare  the  constriction  of  the  normalized  charge  density  profile  for 
50  grid  points  (dashed  curves),  and  for  10  grid  points  (solid  curves). 

It  is  clear  from  the  figure  that  the  spatial  resolution  afforded  by  50 
grid  points  is  adequate  to  produce  accurate  solutions  for  the  constricted 
discharge  at  100  Torr,  while  the  accuracy  of  the  10  point  solution  is 
better  than  expected.  The  constriction  of  the  discharge  in  our  model 
results  from  heating  of  the  neutral  gas  and  the  consequent  outward  flow 
of  neutrals.  This  results  in  a  greater  ionization  rate  at  the  center 
of  the  discharge  than  at  the  wall  due  to  the  exponential  dependence  of 
ionization  on  the  ratio  of  electric  field  to  gas  density.  As  the  charge 
density  increases,  the  discharge  current  also  increases  and  the  axial 
electric  field  must  decrease  if  a  constant  voltage  supply  is  to  be  maintained. 
The  instability  in  Fig.  1  is  thus  limited  by  the  stabilizing  effect 
of  the  external  circuit,  and  a  steadv  state  is  reached.  The  final  approach 
to  steady  state  is  often  accompanied  by  a  slight  expansion  of  the  charge 
density  profile.  The  constriction  is  sometimes  preceded  by  a  semistable 
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plateau  in  the  charge  density  as  occurs  in  Fig.  1  at  about  40  ysecs. 

This  represents  an  initial,  approximate  balance  between  ionization  and 
recombination  before  there  has  been  significant  heating  of  the  gas. 

Results  computed  using  the  spatially  uniform  pressure  model  show, 
the  stabilizing  effect  of  higher  circuit  impedance  and  the  destabilizing 
effect  of  lower  applied  voltage.  We  also  find  that  either  higher  recom¬ 
bination  coefficients  or  higher  pressures  tend  to  accentuate  the  constriction. 
For  pressures  of  1000  Torr,  the  constriction  proceeds  to  such  an  extent 
that  the  spatial  resolution  given  by  50  grid  points  seems  no  longer 
adequate . 

During  the  next  reporting  period  these  techniques  will  be  extended 
to  an  analysis  of  molecular  gases  in  which  attachment  is  important,  e.g., 
the  mixtures  used  in  high  power  CO2  lasers . 


ELECTRON  DENSITY,  n,  electrons/cm 


Figure  1.  Axial  values  of  charge  density,  gas  temperature,  and  gas  density 
for  helium  at  100  Torr  pressure.  Solid  curves  are  solutions 
of  the  complete  problem  for  10  and  50  grid  points.  Dashed  curves 
indicate  solutions  with  spatially  uniform  pressure  and  10  grid 
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Scattering  and  Transport  of  Resonance  Radiation  in  Gases  (Drs.  A.  V.  Phelps 
and  R.  Smick) 

The  objectives  of  this  project  are  to  experimentally  determine  a)  the 
magnitude  and  frequency  distribution  of  radiation  absorbed  and  reemitted 
by  an  atom  for  incident  frequencies  in  the  vicinity  of  a  collision  broadened 
resonance  line,  b)  the  contribution  of  various  scattering  and  excitation 
transfer  processes  to  the  transport  of  the  energy  of  atoms  excited  to  a 
resonance  state,  and  c)  the  range  of  validity  of  various  theories  of  radiative 
and  non-radiative  transport  of  resonance  excitation.  This  project  was 
initiated  in  order  to  improve  the  experimental  data  and  analytical  techniques 
essential  to  the  prediction  of  a)  the  rate  of  radiative  energy  generation 
and  loss  by  plasmas  containing  metal  vapor  and  other  atomic  line  emitters, 
b)  the  intensity  and  frequency  distribution  of  scattered  light,  e.g., 
sunlight  scattered  by  sodium,  atomic  hydrogen,  etc.,  and  c)  the  electrical 
behavior  of  nonequilibrium  discharge  devices  in  which  the  resonance  atom 
density  is  important  in  determining  rates  of  excitation  and  ionization, 
e.g.,  metal  vapor  lasers . 

An  experiment  has  been  assembled  to  measure  the  wavelength  distribution 
of  light  emitted  by  sodium  vapor  which  has  been  illuminated  with  narrowband 
light  from  a  c.  w.  laser.  Sodium  vapor  has  its  lowest-lying  excited  states 
connected  to  the  ground  state  by  optically  allowed  transitions  within  the 
tuning  range  of  c.  w.  dye  lasers  (the  familiar  yellow  D  lines) ,  and  we 
have  chosen  it  as  the  working  substance  in  spite  of  its  tendency  to  degrade 
most  optical  surfaces. 

A  sapphire-windowed  stainless  steel  cell  with  1  mm.  window  separation 
contains  the  sodium  vapor.  In  order  to  accurately  determine  the  vapor 
pressure,  we  use  small  metal  jacketed  thermocouples  and  an  constant- 


temperature  reference.  This  also  houses  stable  preamplifiers  for  the 
monitoring  and  for  feedback  temperature  stabilization  of  the  sodium  cell. 

The  input  laser  power  i.°  normalized  for  changes  in  both  the  pulse 
counting  and  D.C.  outputs  of  the  photomultiplier.  The  dye  laser  will 
produce  output  radiation  narrower  than  .01  nm,  with  output  power  (averaged 
over  10  sec.  internals)  stable  to  +  10%.  The  system  is  presently  being 


calibrated. 


Radiation  Hydrodynamics  (Dr.  D.  G.  Hummer  and  Mr.  P.  Kunasz) 

Line  formations  in  expanding  spherical  geometries.  The  method  discussed 
in  the  last  semi-annual  report  for  solving  the  line  formation  problem  in 
spherical  gaseous  configurations  with  an  arbitrary  radial  velocity  field 
has  now  been  completely  coded,  checked  and  used  to  solve  a  variety  of 
model  problems.  The  results  obtained  by  running  this  code  with  zero  velocity 
are  in  agreement  with  the  results  of  the  variable  Eddington  factor  code 
reported  in  the  last  semi-annual  report.  An  account  of  this  work  has  been 
submitted  for  publication  in  the  Monthly  Notices  of  the  Royal  Astronomical 
Society. 

The  general  conclusion  is  that  velocities  up  to  a  few  times  the  mean 
thermal  velocity  have  a  much  smaller  effect  on  the  source  function  (propor¬ 
tional  to  ratio  of  populations  of  the  upper  to  the  lower  state)  than  does 
the  spherical  geometry.  On  the  other  hand,  even  very  small  velocities 
(1/10  of  the  mean  thermal  velocity)  cause  significant  changes  in  the 
radiation  field.  In  particular,  the  spectral  profile  of  the  radiation  emerg¬ 
ing  from  the  system  is  distorted  in  such  a  way  that  the  long-wavelength 
side  of  the  line  becomes  more  intense,  while  the  short-wavelength  side  is 
reduced  in  intensity.  As  in  static  spherical  models,  the  angular  dependence 
of  the  emergent  radiation  is  very  much  more  drastic  than  for  plane-parallel 
situations . 

As  an  example,  we  consider  an  isothermal  spherical  shell  with  outer 

radius  R  ■  30  times  the  inner  radius  Rc>  The  line  opacity  is  assumed  to 
—2 

fall  off  as  r  and  the  mean  optical  thickness  of  the  shell  is  taken  to 

3 

be  T  =  10  .  A  linear  velocity  law  of  the  form. 

(r-R  ) 

v  (r>  °  y  CR)  UmQ 

is  assumed,  where  v  (R)  is  the  velocity  at  the  outer  radius.  The  source 
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function  in  this  case  is  shown  in  Figure  3,  in  units  such  that  the  Planck 

function  at  line  center  is  unity,  for  maximum  radial  velocities  v  (R)  of  0, 

1  and  2  in  units  of  the  mean  thermal  velocity.  Results  are  shown  for 

-2  -4 

e  i  A  /(C2i  +  A21)  equals  to  10  and  10  *  that  is,  these  systems  are 

very  far  from  local  thfermodynamic  equilibrium.  Here  A2^  and  C2^  are  the 
radiative  and  collisional  deexcitation  rate  coefficients.  For  comparison, 
the  source  functions  for  static,  plane-parallel  models,  R/R^  ■*"  1,  having 
the  same  values  of  e  and  T  are  represented  by  broken  lines.  In  Figure  4, 
the  spectral  profiles  for  these  cases  are  shown.  Here  x  is  the  frequency 
displacement  from  line  center  measured  in  Doppler  units. 

From  these  similar  results  we  can  conclude  that  for  velocities  up  to 
a  few  times  the  mean  thermal  velocity  it  is  satisfactory  for  most  purposes 
to  ignore  the  effects  of  the  velocity  field  in  computing  the  source  function, 
which  is  by  far  the  most  expensive  part  of  the  calculation,  and  to  include 
the  velocity  field  only  in  computing  the  radiation  field  from  the  source 
function. 


•  — 


Log 
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Plasma  Statistics  (Dr.  W.  E.  Brittin) 

Dr.  W.  E.  Brittin  and  A.  Y.  Sahakura  have  discovered  a  general  quantum 

mechanical  transformation  which  transforms  the  usual  description  in  which 

a  system  is  regarded  as  being  made  up  of  particles  (electrons,  atomic  nuclei) 

» 

into  a  description  in  which  the  system  may  be  regarded  made  up  of  bound 
composite  atoms  and  "unbound"  particles.  The  composite  particle  basis 
locks  promising  for  the  treatment  of  both  equilibrium  and  non-equilibriun 
processes  in  many  body  systems.  Both  composite  and  unbound  particles  enter 
the  theory  in  a  rather  simple  manner.  We  expect  to  continue  this  work  and 
apply  it  to  systems  of  interest  in  practice. 


II.  TRANSPORT  COEFFICIENTS  AND  REACTION  RATES 
Program  Leader  -  Dr.  S.  J.  Smith 

The  prediction  of  the  performance  characteri  ics  of  ionized  gas 

devices,  such  as  the  efficiency,  power  output  and  stability  of  high  .power 

* 

lasers,  requires  that  one  be  able  to  calculate  accurate  values  of  elec¬ 
tron  transport  and  reaction  rate  coefficients  for  electrons  and  ions.  In 
order  to  meet  some  of  these  requirements  we  are  engaged  in:  a)  a  determin¬ 
ation  of  sets  of  electron  collision  cross  sections  consistent  with  measured 
electron  transport  and  ionization  coefficients  and  with  recent  electron 
beam  measurements  of  excitation  cross  sections  for  the  upper  excited  states 
of  gases  found  in  molecular  lasers;  b)  measurement  of  cross  sections  for 
molecular  dissociation;  c)  evaluation  and  compilation  of  low  energy  electron 
photon,  and  ion  collision  cross  section  data;  and  d)  determination  of  nega¬ 
tive  ion  characteristics,  such  as  single  and  multiple  p’notodetachnent  cross 
sections,  negative  ion  stability  and  low  electron  energy  scattering  para¬ 


meters  . 
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Electron  Transport  and  Ionization  Coefficients  (Drs.  G.  E.  Chamberlain, 

L.  J.  Kieffer  and  A.  V.  Phelps) 

It  has  been  our  goal  to  compile  sets  of  electron  collision  cross  sections 
consistent  with  electron  transport  and  rate  coefficients  in  order  that  these 

0 

sets  of  data  can  be  used  in  the  modeling  of  weakly  ionized  gas  discharges, 
such  as  high  power  molecular  gas  lasers.  To  this  end  we  have  been  working 
during  the  past  year  on  creating  a  critical  summary  of  the  measurement 
techniques  for  obtaining  electron  transport  and  rate  coefficients  (drift 
velocity,  diffusion,  ionization,  and  excitation). 

We  are  compiling  transport  data  for  ^  in  order  to: 

(a)  compare  the  data  taken  by  differing  techniques, 

(b)  look  for  systematic  effects,  and 

(c)  test  our  criteria  for  selection  of  data  and  try  to  estimate 
the  size  of  systematic  errors. 

Nitrogen  is  being  used  in  these  studies  since  we  believe  most  of  the  problems 
of  measurement,  for  a  non-attaching  gas,  occur  for  and  a  large  body  of 
data  exists  for  N^. 

Currently  our  summaries  of  the  measurement  technique  and  data  collection 
exist  for  drift  velocity  and  transverse  diffusion  (Towns end-Huxley  method). 
Data  has  been  collected  for  ionization,  but  work  on  excitation  has  not 
begun.  Our  recently  completed  review  of  transverse  diffusion  measurements 
emphasizes  the  elementary  need  for  varying  the  geometry,  gas  density,  and 
electron  current  to  check,  for  example,  for  errors  due  to  self  space  charge, 
ion  space  charge,  secondaries  from  ions  or  photons,  and  density  gradients 
at  the  walls.  All  too  often  only  one  point  measurements  are  reported. 

Recent  reports  in  the  literature  have  shown  that  characterization 
of  a  swarm  by  an  anisotropic  electron  diffusion  constant  improves  the 
interpretation  of  high  precision  data  and  removes  the  residual  dependence 
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on  pressure  that  existed  in  some  measurements.  As  a  result  criteria  on 
the  apparatus  geometry  are  available  that  can  be  applied  to  other  earlier 
measurements. 

In  our  current  work  discussions  have  been  put  together  of  the  effect 
» 

of  self  space  charge  on  electron  swarms,  and  a  discussion  of  the  gas  kinetic 
criteria  for  a  swarm  to  be  in  equilibrium.  The  most  universal  criteria 
that  can  be  applied  is  that  the  voltage  V  across  the  drift  space  be  much 
larger  than  the  characteristic  energy  D/p, 

D/p  <<  V. 

Equally  important  is  the  need  to  develop  a  criteria  for  applicability 
of  the  spherical  harmonic  expansion.  This  approximation  is  generally  made 
in  calculating  velocity  distribution  functions  and  relating  cross  sections 
to  transport  coefficients.  Hence  in  calculating  distribution  functions, 
and  in  interpreting  experiments  in  terms  of  the  bulk  parameters,  drift 
velocity,  diffusion,  etc.,  at  high  E/N  (electric  field/gas  number  density) 
it  is  necessary  to  know  the  limits  of  applicability  of  the  spherical 
harmonic  expansion. 

Molecular  Dissociation  Processes  (Dr.  L.  J.  Kieffer) 

No  experimental  work  was  done  under  this  project  during  this  reporting 
period.  Data  taken  last  July  have  been  prepared  for  publication  and  are 
being  submitted  to  Journal  of  Chemical  Physics. 

1.  "Electron  Energy  Dependence  of  the  Energy  and  Angular  Distributions 
of  0  from  Dissociative  Ion  Pair  Formation  in  O^"  R.  J*  Van  Brunt 
and  L.  J.  Kieffer. 

2.  "Breakdown  of  the  Dipole-Born  Approximation  for  Predicting  Angular 
Distributions  of  Dissociation  Fragments"  R.  J.  Van  Brunt. 


and  Mr.  T.  A.  Patterson) 

During  the  present  contract  period  we  have  completed  our  studies 
of  the  photodetachment  of  Se  ,  OH  and  OD  .  The  Se  results  will  be 
published  in  the  Physical  Review;  the  OH  and  OD  studies  will  be  submitted 
to  the  Journal  of  Chemical  Physics.  The  results  of  these  measurements 
have  been  reported  in  previous  Progress  Reports. 

In  recent  months  both  the  tunable  laser  photodetachment  spectroscopy 
device  and  the  fixed  laser  photoelectron  spectroscopy  apparatus  have  been 
utilized  in  studies  of  photodetachment  of  alakali-negative  ions.  All  of 
the  alkalis  have  been  studied  by  at  least  one  of  the  techniques.  The 
resulting  electron  affinities  determined  from  experiment  and  from  theoretical 
techniques  being  developed  by  Dr.  Norcross  at  J1LA  are  summarized  in 
Table  I. 

In  addition  to  these  direct  threshold  measurements,  there  exists  the 

possibility  of  observing  interference  effects  in  photodetachment  at  photon 

energies  corresponding  to  the  opening  of  a  new  channel  which  competes  with 

an  already  open  channel.  This  fact  was  pointed  out  in  detail  by  Norcross 

and  Moores*  in  the  case  of  Na  photodetachment  at  photon  energies  corres- 

2 

ponding  to  the  opening  of  the  lowest  P  exit  channel  in  the  neutral.  From 
a  measurement  of  the  photon  energy  corresponding  to  the  cusp  at  the  channel 

opening  we  can  determine  the  electron  affinities  of  the  alkalis. 

_  2 

Figure  5  shows  Na  photodetachment  data  in  the  region  of  the 

2 

Na(3  Pi/2,3/2)  thresholds.  In  this  case,  the  wavelength  resolution  of  the 
tunable  laser  was  not  adequate  to  resolve  the  individual  fine  structure 
thresholds.  The  cross  section  shape  near  threshold  corresponds  to  a  Wigner 
cusp,  and  is  essentially  in  agreement  with  calculations  performed  by 
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Norcross  and  Moores . J  These  data,  for  various  experimental  reasons,  were 

quite  noisy,  and,  as  a  result,  there  is  no  reason  to  believe  that  the 

oscillations  at  long  wavelengths  are  real.  From  an  analysis  of  these  data 

we  find  EA(Na)  to  be  (0.543  +  .010)  eV. 

The  K  ion  is  a  milch  more  favorable  case  for  experimental  study 

because  the  expected  threshold  falls  in  a  spectral  region  corresponding 

to  a  bright  laser  and  the  spin-orbit  splitting  in  the  excited  state  is 

3 

much  larger  than  for  Na.  Figure  6  shows  K  photodetachment  data  in  the 
2 

vicinity  of  the  ^i/2’3/2  thresholds.  this  case,  the  individual  thresholds 
are  well  resolved  and  we  find  EA(K)  =  (0.5012  +  0.0005)  eV.  That  the 
observed  structure  is  not  a  cusp  is  not  bothersome  since  the  general 
treatment  of  threshold  behavior  predicts  only  a  discontinuity  in  the  derivative 
at  the  channel  opening,  but  not  necessarily  a  derivative  sign  change  at  the 
threshold. 

4 

In  some  very  recent  measurements  of  Cs  photodetachment,  the  threshold 

behavior  is  even  more  unexpected,  as  seen  in  Figure  7.  The  sharp  drop 

o 

in  the  total  Cs  cross  section  near  6682  A  goes,  as  far  as  we  can  tell, 

_3 

to  zero.  Our  present  bound  on  the  minimum  is  3  x  10  of  the  fully  developed 
cross  section,  and  the  shape  of  the  minimum  is  governed  largely  by  the 


far  wings  of  the  laser  line.  Note  that  there  also  appears  to  be  a 

o 

discontinuity  in  the  derivative  of  the  total  cross  section  near  6673  A. 

2 

Similar  behavior  is  seen  near  the  P3/2  threshold,  except  that  the  width 

2 

of  the  hole  is  approximately  a  factor  of  2  greater  than  in  the  case 

and  uhe  minimum  is  not  nearly  so  small.  We  have  no  a  priori  way  of  deciding 


whether  the  channel  opening  occurs  at  6682  X  or  6673  X  (corresponding 


to  electron  affinities  of  0.469  eV  and  0.472  eV  respectively).  Kasdan, 


however,  has  used  the  photoelectron  spectroscopy  apparatus  to  measure  the 


electron  affinity  of  Cs  and  finds  EA(Cs)  =  (0.471  +  0.003)  eV.  His  result 

is  suggestive  that  the  higher  energy  feature  is  the  correct  one  but  not 

conclusive  proof.  Nor  cross**  has  performed  preliminary  three-state  close 

2 

coupling  calculations  on  Cs  photodetachment  in  the  vicinity  of  the  P 


threshold;  these  preliminary  calculations 


that  there  is  a  doubly 


excited  negative  ion  configuration  which  lies  very  near  the  P  threshold. 

Such  a  state  would  give  rise  to  a  Fano  line  shape^  with  a  true  zero  in 

2 

the  total  cross  section.  The  opening  of  the  P  inelastic  channel  would 

interrupt  the  development  of  the  Fano  line  shape  and  give  rise  to  the  higher 

energy  cusp  observed.  If  this  interpretation  is  correct,  then  the  higher 

etergy  cusp  in  Cs  is  understood,  as  is  the  fact  that  the  F3/2  minimum 

is  less  dramatic,  as  well  as  the  qualitative  trend  of  the  threshold  behaviors 

through  the  alkali  ion  series.  As  one  progresses  from  sodium  to  cesium, 

2 

a  temporary  negative  ion  state  appears  and  moves  further  below  the  P 
threshold;  the  development  of  the  resonance  line  shape  corresponding  to 
this  state  is  then  interrupted  by  the  opening  of  the  Inelastic  channel. 

If  this  explanation  is  correct,  then  we  would  expect  to  see  a  Rb  behavior 
which  is  intermediate  between  K  and  Cs  .  Such  an  experiment  is  now  in 
progress.  These  temporary  negative  ion  states  are  frequently  seen  in 
electron  scattering  from  neutrals  just  below  excited  state  thresholds. 

It  may  well  turn  out  that  cesium  represents  a  typical  case,  rather  than 
a  puzzling  exception. 


L 
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TABLE  I. 

ALKALI  ELECTRON  AFFINITIES 

(eV) 

Li 

Na 

K 

Rb 

Cs 

la 

.614 

.537 

.499 

.488 

.466 

lb 

.615 

.540 

.504 

.498 

.481 

2 

.543(10) 

.5012(5) 

.486 (5)p 

•472(5)p 

3 

617 (10)p 

.485(10)p 

.471 (3)p 

4 

^.61 

.53 

^.50 

.45 

^.45 

1  Theory 

-  Norcross : 

a  -  2  state 

calc.,  b  -  3  state  calc. 

(unpublished) 

2  Tunable  Laser  -  Patterson,  et  al  (unpublished) 

3  Photoelectron  Spectroscopy  -  Kasdan,  et  al  (unpublished) 

4  Photodetachment  -  Feldmann,  et  al  (unpublished) 
p  denotes  preliminary  data 
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Fig.  5.  Na  photodetachment  cross  section  in  the  region  of  the 

Na  3  5*1/2, 3/2  thresholds.  The  cusp  occurs  at  the  opening 
of  these  *?  exit  channels.  The  soHd  line  is  a  close 
coupling  calculation  by  Norcross  and  Moores,  convolved 
with  our  experimental  energy  resolution. 
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Fig,  6V  K  photodetachment  data  near  the  2P  thresholds. 
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Cs  photodetachment  cross  section  near  the  2P  thresholds 
The  minimum  near  6682  A  is  broadened  by  our  wavelength 
resolution.  A  possible  interpretation  of  this  resuit  is 
discussed  in  the  text. 
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III.  ENERGY  LOSS  AND  RADIATION  PRODUCTION  BY  ELECTRON  COLLISIONS 

Atomic  and  Radiation  Interaction  Theory  (Dr.  S.  Geltman) 

Free-free  absorption  coefficients.  This  work  is  now  completed  and 
published^.  It  includes  the  calculated  free-free  absorption  coefficients 
for  electrons  and  the  neutral  atoms  He,  C,  N,  0,  Ne,  Ar,  Kr,  and  ;"e,  in 
the  ranges  X  *  .5  to  20  pm  and  T  =  500  -  20000°  K. 

Angular  distributions  for  ionization.  The  mathematical  problems 
encountered  in  applying  the  Coulomb-projected  Bom  approximation  to  the 
ionization  of  atoms  by  electron  impact  have  been  overcome.  Calculations 
for  the  triple-differential  ionization  cross  section  of  hydrogen  have 
been  completed,  and  the  method  is  now  being  applied  to  helium. 

Interaction  of  atoms  with  Intense  laser  fields .  A  theory  has  been 
developed  for  the  absorption  of  ultraintense  radiation  by  atomic  systems. 
It  Is  now  being  applied  to  the  bound-free  of  hydrogen  atoms,  and  will 
later  be  applied  to  the  free-free  absorption  problem.  The  usual  N-photon 
resonant  absorption  of  the  veak  radiation  perturbation  theory  is  replaced 
by  a  non-resonant  non-linear  absorption  in  which  a  monochromatic  radiation 
field  produces  a  continuous  distribution  of  photoionization  electrons. 


^S.  Geltman,  J.  Quant.  Spectrosc.  Radiat.  Transfer  J^3,  601  (1973). 


k. 


-24- 


Theoretical  Atomic  Physics  (Drs .  D.  G.  Hummer  and  D.  W.  Norcross) 

Electron  excitation  cross  sections  for  atmospheric  ions.  The  elaborate 

atomic  structure  models  for  the  Be-like  ions  C  III,  N  IV  and  0  V  were 

used  in  preliminary  testing  runs  of  the  Distorted  Wave  (DW)  collision  pro- 

» 

gram  for  the  purpose  of  establishing  convenient  data  handling  procedures. 

The  present  application,  involving  many  more  reaction  channels  than  in 
earlier  work,  would  involve  the  handling  of  thousands  of  physical  card 
images  if  previous  procedures  were  followed.  Programs  developed  to  pool, 
catalogue,  sort  and  retrieve  the  output  of  the  DW  code  for  use  in  the  program 
SIMMEG,  which  completes  the  calculation  of  excitation  cross  sections,  were 
therefore  tested  and  implemented. 

Sui  sequent  preliminary  runs  of  the  DW  program  for  the  larger  scattering 
problems  were  then  undertaken  and  showed  that  further  alterations  in  the 
program  were  required  due  to  excessive  and  unnecessary  requirements  for 
computer  storage.  These  alterations  are  now  being  made,  in  cooperation 
with  the  Department  of  Physics  of  the  University  of  London,  where  the 
programs  were  originally  developed.  Full  scale  production  runs  of  the 
programs  should  begin  with  the  next  month. 
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Electron-Positive  Ion  Recombination  in  High  Power,  CO^,  Lasers  (Dr.  A.  V. 
Phelps) 

As  the  result  of  r;1  r_.it  experimental  evidence  by  Denes  and  Lowke* 

that  a)  the  "commercial  grade"  gases  used  in  high  power  CC^  laser  contain 

a  fraction  of  one  per  cent  of  ^0  and  b)  the  effective  recombination 

_8  3 

coefficient  for  typical  laser  discharge  conditions  is  less  than  10  cm  /sec, 

we  have  developed  a  model  for  the  positive  ion-molecule  reactions  and 

electron-positive  ion  recombination  processes  occuring  in  high  power  C02 

lasers  using  comercial  grade  gases.  This  model  is  consistent  with  the 

measurements  of  electron  loss  by  recombination  in  under  e-beam  laser 

2 

conditions  by  Douglas -Hamilton  and  the  measurements  of  electron-positive 
ion  recombination  for  H^0+  by  Walls  and  Dunn,  as  well  as  with  the  discharge 
conditions  of  Denes  and  Lowke.  The  principal  features  of  this  model  are: 

1)  N^+  and  C02+ions  produced  by  electron  impact  ionization  of  N2  and 

C02  react  with  H20  to  form  H^0+  in  times  of  the  order  of  0.1  ysec  for 

0.5%  H20  in  typical  laser  mixture  at  atmospheric  pressure.  The  N2+sequence 

A  "f 

has  been  observed  by  Good,  Durden  and  Kebarle  .  The  C02  sequence  is 

predicted"*  to  be  fast  but  has  not  been  observed  experimentally.  The 

■f  <4  *4"  *4“ 

clustered  N,  ion  also  leads  to  H_0  .  The  H_0  ions  will  form  clusters 
A  3  3 

with  H20,  with  the  size  of  the  cluster  depending  upon  the  H20  density,  the 
gas  temperature  and  the  mean  energy  of  the  ions  as  they  drift  through  the 
laser  gas  mixture. ^ 

3 

2)  The  recombination  coefficient  measured  at  JILA  by  Walls  and  Dunn 
for  monoenergetic  electrons  and  H^0+  decreases  as  the  square  of  the  electron 
energies  for  energies  above  about  0.1  eV.  These  authors  have  shown  that 

their  electron  beam  data  can  be  made  consistent  with  the  afterglow^  and 

8  9 

shock  tube  experiments  and  with  theoretical  threshold  laws  by  assuming 
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-k 

that  the  recombination  coefficient  varies  as  (energy)  for  electron 

energies  below  about  0.1  eV.  We  have  made  estimates  of  the  effects  of 

departures  of  the  electron  energy  distribution  from  a  Maxwellian  using  the 

distribution  functions  calculated  previously^  for  ^  and,  recently,  for 

*  11 

various  CO2  laser  mixtures  .  We  also  used  the  approximation  that  because 
of  the  rapid  decrease  in  the  recombination  coefficient  with  increasing 
electron  energy,  only  electrons  at  the  very  lowest  electron  energies  are 
able  to  recombine  with  the  H^O*,  i.e.,  that  the  relative  recombination 
coefficients  are  proportional  to  the  relative  number  of  electrons  with 
energies  near  zero  energy.  Figure  8  shows  a  plot  of  these  relative 
recombination  coefficients  as  a  function  of  average  electron  energy  for 
the  calculated  electron  energy  distributions  in  ^  (solid  curve)  and  for 
Maxwellian  energy  distributions  (broken  curve).  Approximate  recombination 
coefficients  for  electron  energies  well  above  0.1  eV  were  calculated  using 
the  Walls  and  Dunn  data  and  can  be  obtained  using  the  right  hand  scale 
of  Figure  8.  Note  the  more  rapid  variation  with  average  electron  energy 
of  the  recombination  coefficients  calculated  using  the  electron  energy  dis¬ 
tributions  for  N£.  A  somewhat  similar  effect  is  calculated  for  typical 
laser  mixtures  for  the  higher  mean  electron  energies . 

3)  The  electron-positive  ion  recombination  coefficients  measured  by 
Douglas-Hamilton  for  "commercial  grade"  ^  are  shown  by  the  points  of 
Figure  8.  These  data  show  very  nearly  the  same  variation  with  electron 
energy  as  does  the  solid  curve  for  electrons  in  N^.  However,  the  experimental 
values  are  approximately  a  factor  of  3  above  the  values  predicted  for 
H^0+  and  electrons  in  N^.  This  comparison  suggests  the  possibility  that 
the  positive  ions  in  Douglas-Hamilton's  experiment  are  anc* 

that  the  electron  energy  dependence  of  their  recombination  coefficients 


is  similar  to  that  for  This  proposal  requires  that  the  magnitudes 

of  the  recombination  coefficients  for  the  clustered  H^O"*"  ions  be  consistently 

larger  than,  for  the  unclustered  ion,  as  is  observed  for  low  energy  electrons^. 

It  should  be  noted  that  as  the  electric  field  to  gas  density  ratio  E/N 

and  the  electron  energ^  are  increased  the  number  of  water  vapor  molecules 

4*  6 

clustered  about  the  U^O  ion  will  decrease  .  When  this  occurs  we  expect 
the  recombination  coefficient  to  approach  that  of  the  solid  curve.  Measure¬ 
ments  of  ion  breakup  in  the  helium  rich  laser  gas  mixtures  are  not  available. 

4)  If  the  approximate  recombination  coefficients  given  in  Figure  8 

are  extrapolated  to  the  mean  energies  appropriate  to  the  discharge  experiments 

in  of  Denes  and  Lowke,  then  one  predict  a  recombination  coefficient 
+  —8  3 

for  H^O  ions  of  about  2  X  10  cm  /sec  compared  with  a  value  equal  to  or 
_8  3 

less  than  10  cm  /sec  obtained  from  epxeriment.  A  slightly  smaller  value 
is  predicted  for  H^o"*"  and  electrons  in  typical  laser  mixtures,  e.g., 

1:2:3  of  CO^^^^e,  at  the  discharge  E/N. 

5)  A  potentially  important  effect  which  we  have  neglected  in  Figure  8 
is  the  reduction  in  the  effective  recombination  coefficient  resulting 
from  the  depletion  of  the  electron  energy  distribution  function  at  low 
energies  because  of  the  preferential  loss  of  very  low  energy  electrons 

by  electron-H^O^  recombination.  From  the  results  of  Walls  and  Dunn  and 

our  electron  energy  relaxation  calculations^  we  expect  this  effect  to 

become  significant  in  discharges  when  the  fractional  ionization  exceeds 

about  2  X  10  \  This  fractional  ionization  is  about  an  order  of  magnitude 

higher  than  the  densities  reached  in  the  measurements  for  N^  by  Denes  and 

Lowke,  so  that  the  reduction  in  the  recombination  coefficient  is  probably 

small.  Since  the  frequency  of  energy  exchange  collisions  in  mixtures 

12 

containing  significant  amounts  of  CO^  is  much  larger  than  for  N^,  an 


even  larger  fractional  ionization  would  be  required  for  most  CO^  laser 

mixtures  in  order  reduce  the  recombination  coefficient  appreciably.  More 

accurate  estimates  of  this  effect  will  require  a  more  accurate  measurement 

of  the  low-energy  portion  of  the  electron-H,0+  or  H,0+(H„0)  recombination 

»  j  z  n 

coefficient.  * 

6)  The  significant  aspects  ot  i.he  low  electron  recombination  coefficients 
found  by  Denes  and  Lowke  and,  hopefully,  explained  by  our  model  are  that, 

a)  one  understands  the  near  independence  of  the  operating  field 
strength  of  the  glow  discharge  phase  of  laser  discharges  on  discharge 
current  and  can  make  better  predictions  of  the  discharge  excitation  efficiency, 
impedance  characteristics,  etc. 

b)  the  unexpectedly  low  values  of  the  electron-positive  ion  re¬ 
combination  coefficient  for  H^0+  at  electron  energies  characteristic  of 
discharges  suggest  the  possibility  that  other  ions  can  be  found  which  will 
have  more  desirable  recombination  coefficients  than  those  found  in  the 
pure  gases.  Thus,  a  reduction  in  the  electron-positive  ion  recombination 
coefficient  in  a  laser  in  which  the  ionization  is  maintained  by  an  external 
electron  beam  or  by  an  ultraviolet  photon  source  results  in  a  corresponding 
reduction  in  the  energy  required  to  maintain  the  ionization.  It  is  suggested 
that  hydrogen  rich  ions,  as  in  the  NH^  and  H^o"*"  studied  by  Walls  and  Dunn, 
are  more  likely  to  show  unusual  energy  dependences  for  the  electron-positive 
ion  recombination  coefficients  because  of  the  relative  scarcity  of  potential 
energy  curve  crossings  for  the  neutral  and  ionized  species. 

Obviously,  all  of  this  model  is  subject  to  revision  when  the  crutial 
mass  spectrometric  studies  of  ions  are  carried  out  in  high  pressure ,  gas 
laser  mixtures  and  when  more  data  becomes  available  as  to  the  purity  of 
gases  present  in  laser  discharge  experiments. 


RELATIVE  NUMBER  OF  ELECTRONS  AT  ZERO  ENERGY 


RECOMBINATION  COEFFICIENT,  cnT/sec 
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IV.  GENERATION  AND  ABSORPTION  OF  RADIATION 
Program  Leader  -  Dr.  A.  C.  Gallagher 

The  objectives  of  this  program  are  to  determine  some  of  the 
characteristics  of  alkali  metals  and  alkali  metal-rare  gas  systems  of 
importance  in  connection  with  recently  proposed  tunable  lasers*  utilizing 
the  spectra  produced  as  the  result  of  radiative  transitions  in  molecules 
having  a  dissociating  ground  state.  The  projects  being  carried  out  include 
the  determination  of  the  intermolecular  potentials  for  various  alkali 
metal-rare  gas  molecules  and  the  cross  sections  for  electron  impact 
excitation  of  the  resonance  line  of  alkali  metal  atoms. 


*A.  V.  Phelps, 


"Tunable  Gas  Lasers  Utilizing  Ground  State  Dissociation, 


JILA  Report  110.  September  1972. 


A.  C.  Gallagher  and  C.  G.  Carrington) 

The  objectives  of  this  project  are  to  develop  and  use  techniques  for 
the  determination  of  interaction  potentials  for  excited  states  of  diatomic 
molecules  formed  from  resonance  states  and  for  which  the  ground  state  is 
repulsive  or  only  slightly  attractive.  Molecular  states  of  this  type 
determine  the  spectral  absorption  and  emission  profiles  and  the  rates  of 
quenching  of  excited  atoms  and  of  depolarization  of  scattered  light  when 
metal  vapors  are  used  in  a  high  density  of  a  buffer  gas.  An  exciting  appli¬ 
cation  of  the  radiative  properties  of  such  states  is  the  far  UV  laser 
reported  by  Basov  and  coworkers.  Our  recent  measurements  of  the  spectral 
distribution  of  fluorescence  in  the  extreme  wings  of  the  cesium  and  rubidium 
resonance  lines  in  the  presence  of  inert  gas  at  high  densities  have  been 
used  to  obtain  the  interaction  potentials  for  the  ground  state  and  the 
first  two  excited  states  of  the  cesium-  and  rubidium- inert  gas  molecules. 

The  extreme  wing  fluorescence  spectra  is  analyzed  using  the  quasi-static 
model  of  line  broadening  and  readily  yields  much  more  detailed  ini'>rmation 
than  the  previously  available  measurements  of  half-width  and  line  shift 
obtained  near  the  line  center.  In  addition  to  data  yielding  interaction 
potentials,  measurements  of  the  variation  of  the  extreme  wing  intensities 
at  low  inert  gas  densities  contain  information  as  to  the  process  of  formation 
of  excited  molecules  in  bound  states  as  the  result  of  three-body  collisions. 

During  this  reporting  period  we  have  developed  and  tested  a  new  method 
for  learning  molecular  potentials  from  continuum  spectra.  This  method 
is  particularly  useful  for  ground-state  dissociating  molecules  such  as 
Xe^  and  LiXe.  The  molecular  potentials  and  radiative  properties  of  these 


molecules,  which  are  obtained  by  this  technique,  are  major  determinants 
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of  their  behavior  in  molecular-dissociation  laser  systems. 

The  new  method  for  unique  determination  of  the  potential  compliments 
the  method  previously  developed  under  this  contract.  Whereas  the  previous 
method  utilized  the  temperature  dependence  of  high-pressure  emission  spectra, 

9 

the  new  method  uses  the  pressure  dependence  at  one  temperature.  Thus  it 

is  particularly  useful  for  materials  such  as  LiXe,  where  the  temperature- 

dependence  method  requires  very  high  temperatures  with  concommitant 

problems  of  containing  the  highly-reactive  lithium.  A  detailed  description 

of  the  new  method  is  being  prepared  for  publication. 

The  obvious  advantage  of  a  second  potential-determination  method 

is  the  opportunity  to  check,  the  self-consistency  of  potentials  obtained 

by  two  methods.  We  have  made  such  a  check  for  the  RbXe  ground-state 

2  2 

dissociating  molecule.  Each  method  obtains  a  set  of  A  and  X  potentials; 
the  two  sets  are  compared  in  Figure  9.  The  differences  between  the  two 
sets  of  potentials  are  within  the  combined  experimental  errors  of  the  two 
experiments  and  both  sets  of  potentials  yield  the  same  high-pressure 
RbXe  spectrum  at  300°  K.  The  reasonable  agreement  between  these  two  sets 
of  potentials  gives  us  confidence  in  the  models,  while  the  remaining 
differences  point  out  the  need  for  refinements  in  the  experimental  technique. 


Electronic 
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Excitation  of  Metal  Atom  Resonance  Lines  (Dr.  A.  C.  Gallagher  and  Mr. 

D.  Leep) 

Measurements  of  the  cross  section  for  electron  excitation  of  the 
first  excited  state  of  Li  have  been  completed  and  are  being  prepared  for 
publication.  Figure  ltf  shows  a  comparison  of  these  results  with  one  of 
the  theoretical  calculations^-.  More  recent  theoretical  calculations 
yield  even  larger  cross  sections  near  the  maximum.  Experimental  data 
such  as  this  is  necessary  for  the  prediction  of  ele<  tror  excitation  rates 
in  the  proposed  LiXe  laser  as  well  as  for  testing  the  close  coupling 
technique  being  used  for  the  calculation  of  excitation  cross  sections 
for  large  numbers  elements  of  light  mass  and  high  chemical  reactivity, 
e.g.,  atomic  oxygen.  These  results  show  the  continued  importance  of 
experimental  measurements  of  electron  scattering  from  neutral  atoms  of 
technological  significance. 
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Figure  10.  Electron  excitation  cross  sections  for  excitation  of  the  2s 
state  of  lithium  as  measured  using  the  670.7  ym  resonance 
line  (solid  curve)  and  as  calculated  by  Burke  and  Taylor 
(dashed  curve) . 
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Plasma  Spectroscopy  (Dr.  J.  Cooper  and  Mr.  J.  Baur) 

Work  with  the  pressure  driven  shock  tube  (J.  Baur)  with  cesium  alum 
added  as  a  powered  aersol  indicates  that,  at  a  temperature  of  6000°  K, 
electron  densities  of  up  to  5  X  10  cm  can  be  easily  obtained.  Measure— 
ments  of  electron  density  have  been  performed  using  a  two  wavelength  inter¬ 
ferometer  (in  a  Michelson  arrangement)  and  also  using  the  Stark  broadening 
of  H-o.  Consistency  of  the  measurements  are  to  about  10%.  The  Stark 
b) oadening  of  both  the  calcium  II  H  &  K  lines  and  the  sodium  D  lines  are 
being  measured.  At  the  low  temperatures  considered  here,  the  Stark  broaden¬ 
ing  is  increasing  rapidly  with  decreasing  temperature.  For  calcium  II 
o 

3968  A,  the  ratio  of  width  to  electron  density  at  5000°  K  is  roughly  twice 
that  at  20,000°  K,  which  is  consistent  with  theoretical  predictions. 

Theoretical  work  has  been  performed  on  several  aspects  of  resonance 
broadening.  It  Is  found  that  coherent  transfer  of  excitation  from  one 
excited  atom  to  another  ground  state  atom  leads  to  a  shift  of  the  resonance 
3-Ine  itself,  but  for  radiation  terminating  on  the  upper  level  of  the  resonance 
transition  there  are  no  additional  line  broadening  effects  (ocher  than  those 
included  in  the  traditional  impact  theory).  When  two  identical  atoms  (a 
fixed  distance  apart)  interact  via  their  radiation  field  it  is  possible  to 
form  both  a  symmetric  and  an  antisymmetric  combination  of  states.  At  small 
distances  the  symmetric  (super  radiant)  state  decays  with  twice  the  rate 
associated  with  an  isolated  atom,  whereas  the  decay  of  the  antisymmetric 
state  is  very  small.  We  have  examined  the  role  of  relative  motion  (Doppler 
effect)  on  this  coherent  interaction,  and  find  that  the  radiative  coupling 
Is  reduced  to  essentially  zero  so  that  these  states  decay  as  for  isolated 
atoms  when  the  atoms  can  move  a  distance  of  roughly  \/2  in  their  natural 


lifetime. 
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